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ABSTRACT 

The center-to-limb variation (CLV) describes the brightness of the stellar disk as a function of the limb angle. Across strong absorption 
lines, the CLV can vary quite significantly. We obtained a densely sampled time series of high-resolution transit spectra of the active 
planet host star HD 189733 with UVES. Using the passing planetary disk of the hot Jupiter HD 189733 b as a probe, we study the CLV 
in the wings of the Can H and K and Nai D, and D 2 Fraunhofer lines, which are not strongly affected by activity-induced variability. 
In agreement with model predictions, our analysis shows that the wings of the studied Fraunhofer lines are limb brightened with 
respect to the (quasi-)continuum. The strength of the CLV-induced effect can be on the same order as signals found for hot Jupiter 
atmospheres. Therefore, a careful treatment of the wavelength dependence of the stellar CLV in strong absorption lines is highly 
relevant in the interpretation of planetary transit spectroscopy. 
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1. Introduction 

The discovery of the first transiting extrasolar planet 
HD 209458 b ( |Charbonneau et ar]|2000| l followed five years 
after the detection of the first extrasolar planet around a main- 
sequence star (51 Pegb, Mayor & Queloz|1995| l. Both 51 Pegb 
and HD 209458 b are Jovian planets orbiting their host stars 
at a distance of only a few hundredth of an AU, so-called hot 
Jupiters. At these distances, the atmospheres of the planets are 
exposed to strong stellar radiation fields. Unequaled in the solar 
system, the atmospheres of hot Jupiters can best be studied in 
transiting exoplanetary systems, where they can be observed in 
absorption against the bright background of the stellar disk. 

In particular the stellar extreme ultraviolet and X-ray ra¬ 
diation efficiently deposit energy in the upper planetary atmo¬ 
sphere, where it is thought to derive a Parker-like planetary wind 
(e.g., [Watson et al.| 1981| l. The wind material reaches beyond 
the optical radius of the planet and effectively enlarges the ra¬ 
dius of the disk, especially at wavelengths were the absorption 
and scattering cross section of the material is large, such as the 
hydrogen Lyman a line. Expanding planetary atmospheres like 
this have, indeed, been detected by means of ultraviolet transit 
spectroscopy, for instance, in the hot Jupiters HD 189733 b and 
HD 209458 b ( [Vidal-Madjar et al.|2003[ Lecavelier Des Etangs 
|et al.|2010| ). Eurther ultraviolet spectral features attributed to the 
upper planetary atmosphere have been found in HD 209458 b 
and WASP-12 b (e.g., Linsky et al. 2010[ Eossati et al. 2010 1 . 
Even an X-ray detection of the atmosphere of HD 189733 b has 
been reported ( [Poppenhaeger et al.|2013| l. 

The planetary atmosphere can also be probed by in-transit 


excess absorption in the optical regime (e.g., Seager & Sasselov 


|2000[ [Brown|2001] ). Variations in the planetary radius due to the 
presence of an atmosphere can be probed, e.g., in the Eraun- 
hofer Nai Di and D 2 lines at 5900 A. In 2002, Charbonneau 


* Based on observations made with UVES at the ESO VLT Kueyen 
telescope under program 089.D-0701(A). 


|et al.| announced the detection of in-transit excess absorption at 
a level of (2.32 + 0.57) x 10“^ in the Nai Di and D 2 lines of 
HD 209458 observed with the “Space Telescope Imaging Spec¬ 
trograph” (STIS) on board the Hubble Space Telescope (HST). 
While the first detection of atmospheric sodium was achieved 
with the HST, the visual regime is also accessible with ground- 
based instrumentation. 

After its space-based discovery, the atmosphere of 
HD 209458 b was also the first to be detected with ground-based 
instrumentation. In a careful reanalysis of a data set earlier 
discussed by [Narita et al.| P005| l, [Snellen et al.| ( [2008 1 found 
excess absorption in the Nai Di and D 2 lines of HD 209458b 
at a level of (1.35 + 0.17) x 10“^ in 0.75 A wide spectral 
bands centered on the Nai Dj and D 2 line cores. The signal is 
less pronounced when broader spectral bands are used. Today, 
ground-based transmission spectroscopy is regularly used to 
study planetary atmospheres; a number of successful detections 
of atmospheric sodium have been reported in the literature (see, 
e.g., [Snellen et alT 2008[ Wood et al. 201 1[ [Sing et aL][2012[ 
[Zhou & Bayliss|2012[ Burton et al.|2^^ 

The hot Jupiter HD 189733 b is one of the best-studied plan¬ 
ets beyond the solar system to date. The planet orbits a K-type 
host star at a distance of 0.03 AU and transits the stellar disk 
every 2.2 d ( [Bouchy et al.[[2005| ). In contrast to HD 209458, 
HD 189733 is a highly active star. The spectrum shows promi¬ 
nent chromospheric cores in the Ca ii H and K lines and, conse¬ 
quently, a large S-index of 0.525 ([Wright et al.[[2004| [Knutson] 
et al.|2010[ l. This is consistent with rotational photometric vari¬ 


ability on the level of a few percent? and an X-ray luminosity 


of 1.3 X 10^® ergs ‘, which makes HD 189733 one of the most 
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X-ray luminous K-type stars in the solar neighborhood, and cer¬ 
tainly puts it among the most active planet host stars known to 


date dSchmitt et al.|1995[ Croll et al.|2007[ Poppenhaeger et al. 
|2013| l. Its distance of only 19.3 pc and resulting apparent bright¬ 

ness of y = 7.67 mag have made it a prime target for detailed 
follow-up campaigns at all wavelengths ranges. 

HD 189733 was the target of an HST campaign carried out 
with the Advanced Camera for Surveys (ACS) presented by Pont 
et al.|(|2008|l. The authors find an “al most featureless” spec trum 

(2012 1 pre- 


Huitson et al. 


between 5500 and 10500 A. Later, 
sented another set of transmission spectra of HD 189733 b oh 
tained with HST/STIS. In their analysis, [Huitson et al.| ( |20T^ 
found excess absorption at a level of (9 + l)x 10“"^ inthe Na i Di 
and D 2 lines. The excess absorption is concentrated in the line 
cores, which [Huitson et al.| ( |2012| l attribute to planetary high- 
altitude haze or a substantially subsolar sodium abundance of the 
planet. Both Pont et al. ( |2008| l and Huitson et al. ( |2012| l found 
clear evidence for starspot occultations in their light curves, 
which underlines the role of the high level of stellar activity in 
HD 189733. 

[Redfield et al.| ( |2008| ) reported on Nat excess absorption in 
HD 189733, which is about twice as strong as that observed in 
HD 209458 b by [Snellen et al.| ( (2008] ). |Redfield et al.| ( p008l l used 
a different convention to quantify the excess absorption; a con¬ 
version has been carried ou t by [Snellen et ah] (|2008|l. I n a re¬ 
analysis of the data used by Redfield et al. Jensen et al. ( |2011 1 
confirmed their results and reported excess absorption at a level 
of (5.26 + 1.69) X 10“"^ across a 12 A wide spectral band. 

Recently, I Wyttenbach et al.|p015| l reported on another mea¬ 
surement of excess absorption at a level of (3.2 + 0.3) x 10“^ 
in the Nat Di and D 2 lines of HD 189733b based on transit 
spectroscopy obtained with HARPS. The quoted number refers 
to the same 0.75 A wide bands used earlier by Snellen et al. 
( |2008| l. [Wyttenbach et al.| ( |2015[ ) conclude that the major part of 
their signal is produced in the line cores and could even quan¬ 
tify the width of the signal, giving a full width at half maxi¬ 
mum (FWHM) of (0.52 + 0.08) A. Additionally, there seems to 
be slight blue shift, which may be indicative of a wind in the 
planetary atmosphere of HD 189733 b blowing from the day to 
the night side, and is similar to what is observed in the case of 
HD 209458 b ( [Snellen et al.TIOTOl l. 

While transit spectroscopy is a powerful technique for study¬ 
ing the planetary atmosphere, it is only seen in absorption against 
the stellar disk, whose appearance is, therefore, equally relevant 
for an appropriate interpretation of the results. A wavelength- 
dependent change in the brightness of the disk from the center to¬ 
ward the limb has long been observed on the Sun (e.g.. Pierce & 
[Slaughter! 1977[|Neckel & Labs|1994] l. The center-to-limb vari¬ 
ation (CLV) is attributable to a limb-angle dependent viewing 
geometry, effectively exposing ever higher layers of the solar at¬ 
mosphere to the observer targeting the limb. As higher layers are 
usually cooler and emit less light in the optical, the effect is of¬ 
ten more specifically dubbed limb darkening. This term mostly, 
but not always, describes the appearance of the stellar disk in 
the visual band. At other wavelengths, such as in the X-ray and 
ultraviolet regime, the solar atmosphere can show limb brighten¬ 
ing (e.g., Schlawin et al. [20i0{ [Poppenhaeger et al.|2013)[Llama[ 


|& Shkolnik[M5F 

In the analysis of transit light curves of extrasolar planets, 
the CLV has long been acknowledged as a critical factor (e.g., 
Muller et al. 2013[l. In the modeling, the brightness distribu¬ 


tion of the stellar disk is usually parametrized by one of many 
available so-called limb-darkening laws. Among the most fre- 
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quently applied parameterizations are the linear and quadratic 
limb-darkening laws, according to which the limb-angle depen¬ 
dent intensity is given by 


7linear(/^) 

= 1 ) 


1 - a (1 - p) 


and 


7quadratic(/^) 


/(/r= 1) 


l-a(l-p)-b(l-pf , 


( 1 ) 


( 2 ) 


where a and b are the linear and quadratic limb-darkening 
coefficients, and p is the limb angle (see Sect. 113 . Limb- 
darkening coefficients for a wide range of stellar parameters, 
limb-darkening laws, and photometric bands have been tabulated 
(e.g., [Claret|2^ . 

While the tabulated coefficients usually refer to broad photo¬ 
metric bands, the wavelength dependence of the CLV in the stel¬ 
lar atmosphere can be very strong on shorter scales. In the opti¬ 
cal, the wavelength dependence is strongest across the profiles of 
deep stellar absorption lines, such as the Caii H and K or Nai Di 
and D 2 lines (e.g., Athay et al.[1972[ l. Yan et al.] ( 2015[ l presented 
a study on the CLV in the solar Fraunhofer lines observed in 
solar eclipse spectra reflected by the Moon. The authors clearly 
detected the CLV of the solar atmosphere in their spectra and 
argue that the effect becomes most prominent if narrow spectral 
bands covering the strong (Fraunhofer) lines are used. 

In fact, strong modulations in the stellar CLV over small 
wavelength ranges also dubbed differential CLV (or differential 
limb darkening) are a known problem in the analysis of the Nat 
features in extrasolar planets. For instance, [Charbonneau et ah] 
P002[ l provide a careful discussion on the effect to rule out that 
their signal is caused by the CLV. Also Sing et al. ( 2008| l discuss 
the CLV in HD 209458 extensively and present CLV-corrected 
transit light curves. Among others, the authors conclude that 
limb darkening in the stellar lines is smaller than in the con¬ 
tinuum and that one-dimensional stellar atmospheres can repro¬ 
duce the observed effects. In fact, many authors, including some 
of those referenced earlier, have discussed the issue so that the 
above list is by no means complete. 

The treatment of the CLV is further complicated by the avail¬ 


ability (or lack thereof) of different model predictions. Claret 
( |2004[ l give both limb-darkening coefficients derived from (one¬ 
dimensional) PHOENIX and ATLAS model atmospheres, which 
are similar but not identical. [Muller et al.[ ([201 3|l find both 


sets of coefficients to be consistent with Kepler data. Hayek 
[et al.[ ( |20l3 compare CLV predictions resulting from one- and 
three-dimensional model atmospheres to transit light curves of 
HD 189733 and HD 209458. For the latter, they find a consider¬ 
ably better match with the three-dimensional model. While this 
is also compatible with their finding for HD 189733, their data 
remained insufficient for a conclusive result in this case. 


In this work, we present the analysis of the CLV in the Ca ii 
H and K and Nat Di and D 2 lines in HD 189733 and discuss its 
relevance for searches of atmospheric excess absorption in the 
Nai lines via transit spectroscopy. In Sects. [^to|2 we introduce 
the formalism and models used in our analysis. Our observa¬ 
tions and the data analysis are presented in Sects. [^to[^ Finally, 
we discuss the results in Sect. and present our conclusions in 
Sect.il 
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2. Synthetic data 

Throughout this work, we present synthetic spectra and simu¬ 
lated light curves of transiting planetary systems. Here, we out¬ 
line the simulations and input data. 


3. Transit iight curves and their combination 

Transit light curves, whether derived from broadband photome¬ 
try or high-resolution spectroscopy, are among the most impor¬ 
tant tools to study exoplanets. Below, we briefly discuss the most 
relevant aspects required in our work. 


2.1. Limb-angle dependent spectra 

The limb angle, 6, is the angle between the outward normal of 
the stellar atmosphere and the vector pointing from the center of 
the star toward the observer. Usually its cosine p - cos(0) is used 
to specify the limb angle. With this convention, p - 0 refers to 
the limb and yu = 1 to the center of the stellar disk. 

All synthetic spectra used in this work are based on Kurucz 
model atmospheres with solar metallicity ( [Castelli & Kuruc^ 
2004[ Kurucz|1970|l. The spectra themselves were generated us 


ing the spectrum prograrrj^by R.O. Gray (Gray & Corbally 


|1994| l, which facilitates the generation of both disk-integrated 
stellar spectra and limb-angle dependent specific intensities. 
Based on the assumption of plane-parallel atmospheres and local 
thermodynamic equilibrium (LTE), spectrum first computes the 
number densities of electrons and all relevant species in the at¬ 
mosphere and solves the equation of radiative transfer to obtain 
a synthetic spectrum. Information on individual spectral lines, 
viz., their wavelength, the ion, the lower and upper energy lev¬ 
els, and the product of statistical weight and oscillator strength 
are specified via a line list file. In our calculations, we relied on 
the line list shipped along with the program. 

Whenever we use limb-angle dependent specific intensities, 
we generate 20 spectra with limb angles between p - 0 and 
p - I, i.e., a spacing of Ap = 0.05. To avoid numerical problems 
at the very edge of the stellar disk {p - 0), we replace the lower 
limit hy p - 0.001 in the calculations and assume that the spec¬ 
trum completely vanishes for p -0. The impact of this approxi¬ 
mation remains small because the stellar atmosphere is seen at a 
very steep angle so close to the limb and only contributes a tiny 
fraction to the total observed flux. Spectra at intermediate limb 
angles are linearly interpolated based on the spectra on the grid. 
All spectra are equidistantly sampled on an equidistant wave¬ 
length grid with a spacing of 0.01 A per bin. 


2.2. Discretized stellar surface 

To simulate the spectra of a stellar disk, potentially partially 
occulted by a planet, we use a discretized stellar surface with 
elements defined according to the prescription by |Vogt et 3T] 
( ]1987| l. In particular, we use 1001 latitudinal rings and a total 
of 250000 surface elements. In the Vogt et al. discretization, 
the surface elements have approximately the same size, which 
is about 5 x 10“^ rad or 0.17 square degree in our case. 

The spectrum is obtained by summing the contributions of all 
visible surface elements, assigning to each a spectrum according 
to its projected area, limb angle, and radial velocity. A transiting 
planet is implemented by occulting all stellar surface elements 
covered by the planetary disk. This is equivalent to setting their 
projected area to zero during occultation. 

Given this configuration, the surface element with the largest 
projected area covers the tt x (5 x 10“^)“' = 62 500th part of 
the visible stellar disk. Therefore, a Jovian planet occulting one 
percent of the disk would be sampled by about 625 stellar surface 
elements, which gives an idea of the statistical uncertainty in the 
simulations. 


* http://www.appstate.edu/~grayro/spectrum/spectrum.html 


3.1. The standard planetary system 

We present a number of model calculations referring to what 
we call our standard planetary system, which are used when¬ 
ever we do not discuss a specific planetary system such as 
HD 189733. The standard planet represents a typical hot Jupiter. 
It has one Jovian radius and orbits a star with solar radius in a 
distance of ten solar radii at an inclination of 90 degrees. The 
orbital period is 2 days. The standard star is a slow rotator with 
iisin(0 = 0.1 kms“* seen at an inclination of 90 degrees (i.e., 
equator-on). Assuming the convention that the primary transit 
occurs at orbital phase zero, the transit lasts from orbital phase 
-0.018 to h- 0.018, and the transit duration is 1.69 h. 


3.2. Naming conventions for light curves 

We obtain the time-dependent flux, in a spectral band cen¬ 
tered on some feature of interest and its equivalent, /r(f,), in the 
reference band, where the f, refer to the times of our hypotheti¬ 
cal measurements. Further, we let fjt and foot refer to all measure¬ 
ments in- and out-of-transit. We dub //(f,) and friti) two light 
curves. In transit analyses, these are often normalized by their 
out-of-transit levels, so that 


nxiti) = 


fxitj) 

fxitoot) 


(3) 


where X represents f or r and the bar indicates the mean value. 
We refer to nx as the normalized light curve. The difference, 
of the normalized light curves, i.e.. 


d{ti) = 


//(6) friti) 


ffit oot) fritoot) 


: rifiti) - ririti) , 


(4) 


traces differential changes in between them. We dub dfi) the dif¬ 
ference curve (DC); this form was used by Charbonneau et ah] 
P002|l in their analysis. To quantify the difference in transit 


depth, we follow Charbonneau et ak] ( 2002| l in defining the in¬ 
transit DC excess (DCE) 


DCE = DC(fit) - DC(foot) . (5) 

We note that another natural choice for studying differential 
changes is the ratio, r(f,), of the normalized light curves 


rit.) = X 

/'■(f') JJiiffd nriti) 


( 6 ) 


where we subtract one to obtain the same out-of-transit level as 
for the DC. Yet, the difference between the DC and the ratio 
curve 


DC(f)-r(f) = DC(f)(l-(7) 

is small for typical planetary transit light curves with transit 
depths of a few percent. Therefore, we limit the following dis¬ 
cussion to the DC. 
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3.3. The appearance of the difference curve 

Assume we are given the normalized light curves of a feature, 
LCf(ti), and a number of normalized reference light curves. 
For instance, these light curves could have been ex¬ 
tracted from a series of spectra by integrating the signal in var¬ 
ious bands. In this particular case, we assume that LC/(f,) and 
LCrj „(f,) are normalized transit light curves of our “standard 
planetary system” (s ee Sect.|3.1|l, which we obtain from conven¬ 
tional transit models ( |Mandel & Agol|2002| l. The only difference 
between the feature and reference light curves is the stellar limb 
darkening. For the sake of simplicity and without loss of gen¬ 
erality, we assume the linear limb-darkening law in this demon¬ 
stration (see Eq.[T]l. 

For the feature light curve, LC we specifically assume 
a value of 0.5 for the linear limb-darkening coefficient, a. Then, 
we generate two reference light curves; the first, LCrfi.siti), using 
a limb-darkening coefficient of 0.3 and another one, LCrpniti), 
assuming a value of 0.7. The resulting model light curves are 
shown in Fig. [T] Compared to LCr.o.sCh), the feature shows 
limb darkening, but it shows limb brightening with respect to 
LCrpniti). We emphasize that all of the light curves show limb 
darkening, and limb brightening is only seen as a relative effect. 

In Fig. ^ we show both the DCs obtained by combining the 
feature light curve, LC/(f,), with either of the reference light 
curves. The DCs show a distinct behavior, which differs sub¬ 
stantially from that of the input transit light curves (cf. Fig. [^. 
In the case of identical limb darkening for the feature and refer¬ 
ence light curve, the DC becomes zero. 



Fig. 2. Difference curve, DC(f,), for a limb-darkened (solid, LC/(f,) vs. 
LCrfi.iiti)) and a limb-brightened (dashed, LC/(f,) vs. LCr,oj(ti)) feature. 


When the feature is limb brightened compared to the refer¬ 
ence light curve, the DC shows a drop below the out-of-transit 
level during in- and egress. At the limb of the stellar disk, the 
planet blocks more light from the feature than it blocks from the 
reference. In the middle of the transit, with the planet in the cen¬ 
ter of the stellar disk, the DC shows a peak because the situation 
is reversed. Compared to the feature, more of the reference light 
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is emitted in the center of the stellar disk, where it is now oc¬ 
culted by the planetary disk. When the feature is limb darkened, 
the DC is inverted, but otherwise keeps its main characteristics. 


4. Differential light curves and the stellar disk 

The presence of excess absorption attributable to an exoplane¬ 
tary atmosphere during transit is often derived by analyzing dif¬ 
ferential light curves, such as the DC introduced in Sect. |3.2| 
Usually, the spectral flux is integrated in a band centered on the 
specific feature under consideration and compared to the flux ob¬ 
served in various reference bands (e.g., |Charbonneau et al. 2002 [ 
[Snellen et al.||2008] l. Depending on the characteristics of these 
bands, the limb-angle dependence of the stellar spectrum, i.e., 
the stellar atmosphere, can imprint a strong signal on the DC 
totally unrelated to the planetary atmosphere. 


4. 1. Limb-angle dependence of spectral lines 


X 
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Fig. 3. Normalized stellar spectrum around the sodium doublet for two 
limb angles p = 0.1 (dashed) and p = 1.0 (solid) and three effective 
temperatures. 


In Fig.|^ we demonstrate the effective temperature and limb- 
angle dependence of the Nai D] and D 2 lines using stars with 
effective temperatures of 4000 K, 5000 K, and 6000 K. Irrespec¬ 
tive of limb angle, the width of the Na i Di and D 2 lines increases 
when the effective temperature decreases. This behavior can pri¬ 
marily be attributed to pressure broadening, viz., to the temper¬ 
ature and pressure dependence of the line damping constants. In 
solar-like stars, pressure broadening of the Nai Di and D 2 lines 
is dominated by van der Waals forces (Gray 2008| Fig. 11.4), 
whose influence is stronger in the denser atmospheres of cooler 
stars. 
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At all temperatures, there is a strong dependence of the line 
profile on the limb angle. In the center of the stellar disk, the 
line wings are more pronounced than at the limb. At the stellar 
limb, higher layers with less pressure are seen, where pressure 
broadening is weaker. Although the limb-angle dependence is 
more clearly seen in cooler stars, this is mainly attributable to 
the larger line width. If the wavelength was scaled by the Nai 
line width, the behavior of the line wings would actually become 
quite similar. In contrast to the wings, the depth of the line core 
remains virtually unchanged, which has an immediate and im¬ 
portant consequence: compared to the core and the continuum, 
the line wings of the Nat doublet show limb brightening. In fact, 
not only the Nat lines, but all spectral lines show distinct CLV. 

In Fig. 1^ we show the ratio, RiA), of the specific intensity in 
the center (//=!) and at the limb (ju = 0.05) of a stellar disk with 
an effective temperature of 5000 K and log(g') = 4.5. In terms of 
limb darkening (see Sect. [1]), this ratio corresponds to 


R(A)^ 


7(0.05) 

/(l.O) 


( 8 ) 


While the (quasi-)continuum sections of the spectrum show a 
mostly constant ratio, the spectral lines show distinct behavior. 
Ignoring the sodium doublet for the moment, we find that R{A) 
increases for the majority of spectral lines covered in Fig. ^ 
These lines, therefore, show limb brightening with respect to the 
continuum, however, a number of lines also show the reverse ef¬ 
fect. 

Whether a spectral line shows limb brightening or darkening 
with respect to the continuum depends on its height of forma¬ 
tion in the stellar atmosphere. Lines formed in the upper lay¬ 
ers become stronger on the limb, where their depth and equiva¬ 
lent width increases because the line of sight resulting from the 
steeper viewing angle favors higher layers. Spectral lines formed 
in deeper layers become stronger in the center of the disk, where 
we look deeper into the atmosphere. While the distribution of 
spectral line formation in the stellar atmosphere depends on the 
details of the stratification, it can be said that the height of for¬ 
mation for lines originating from a single ion tends to be higher 
for lower excitation potentials (e.g., Grossmann-Doerth|1994| l. 

The strongest lines, such as the Na I doublet, are formed over 
a wide range of depths in the atmosphere. Therefore, they show 
the most complex behavior of the spectral lines visible in Fig. 
(top). The line cores show immense optical thickness and are, 
therefore, formed high up in the atmosphere. In fact, at/j.- 0.05 
the intensity in the line cores has almost dropped to zero already 
(Fig. 0 bottom). This corresponds to an effective decrease in 
the stellar radius or an increase in the planet-to-star radius ratio 
in the cores of the Nai lines. While the cores, therefore, show 
strong limb-darkening with respect to the continuum emission, 
the line wings exhibit the opposite behavior. The wings are limb 
brightened with the details depending on the exact spectral band 
chosen. Again, we emphasize that the line wings show only rela¬ 
tive limb brightening compared to the continuum. In the bottom 
panel of Fig. we show the normalized specific intensity of the 
line core, a point in the line wing, and a point in the surrounding 
quasi-continuum. The curves clearly demonstrate that the line 
core shows limb darkening and the line wing shows limb bright¬ 
ening relative to the continuum. 

Figure clearly shows that the assumption of a homoge¬ 
neously limb-darkened stellar disk becomes inappropriate as 
soon as individual spectral lines are considered. The limb-angle 
dependence is most prominent in deep and broad lines, such as 
the Can H and K lines or the Nat doublet lines. The latter are 



0 0.2 0.4 0.6 0.8 1 

F 

Fig. 4. Top: the ratio of stellar spectra (Teff = 5000 K, log(g') = 4.5) at 
limb and center (p = 0.05 and /j = 1.0) shown with a solid line. For 
comparison, the scaled central spectrum (dashed line). Bottom: nor¬ 
malized specific intensity in the p = 0.001 - 1 range in the line core 
(5889.95 A, crosses), the line wing (5889 A, squares), and the quasi 
-continuum (5864.6 A, circles). 


a main focus of the study of exoplanet atmospheres (e.g., Red-| 


field et ar]|2008 |Charbonneau et al.]|2002[ [Snellen et al. 1120081 

Huitson et al.|201^ 


4.2. Simulated Na I transit light curves 

Using the techniques outlined in Sect.|^ we simulated spectral 
time series covering a transit of our standard exoplanetary sys¬ 
tem (see Sect. HU- In the simulations, we used a grid of stellar 
effective temperatures ranging from 4000 K to 7000 K in steps 
of 500 K; log(g') was 4.5 in all simulations. 

We simulated the spectra for 60 equidistant time points be¬ 
tween orbital phases of -0.02 and 0.02, where the transit cen¬ 
ter occurs at phase zero. The temporal cadence of the simulated 
spectra, therefore, is (0.04 x 2 d)/60 = 1.92 min. The numbers 
quoted in the following refer to spectra not broadened by instru¬ 
mental resolution, however, we found its impact remains weak 
as long the considered spectral bands cover a few spectral reso¬ 
lution elements. 


4.2.1. Difference curves with broad feature intervals 


In their analysis of the atmosphere of HD 209458 b, jCharbon- 
neau et aT] (2002 | ) used bands covering the entire sodium-line 
doublet. In particular, they defined narrow (n), medium (m), and 
wide (w) bands for the blue and red reference bands and the cen¬ 
tral band covering the sodium lines. 

jCharbonneau et ah] ( |2002| ) detect excess absorption by the 
planetary atmosphere through a negative excess in what we 
dubbed the DC (see Sect. |3.2| l. As the authors point out, the CLV 
is a potential source of this kind of an excess. However, they ver¬ 
ified that it is of no concern in their analysis because the effect 
remains small (see their Eq. 6). 


We used the spectral bands defined by Charbonneau et al. 
( 2002| l and reproduced their DCE estimates (see Sect. 3.21 for 
our simulations of the standard planetary system. Our results are 
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Table 1. Simulated in-transit excess absorption diagnosed by the DC 
for the narrow (n), medium(m), and wide (w) spectral bands defined by 
ICharbonneau et ar]j2002| l. 


Table 2. Difference curve excesses for narrow reference bands with full 
widths of 0.75 A, 1.5 A, and 3 A centered on the Nai Di and Do lines. 
All DCEs are given in units of 10“^. 


Teff 

[K] 

n 

[10-5] 

m 

[10-5] 

w 

[10-5] 

Teff [K] 

0.75 A 

D 2 Di 

1.5 A 

D 2 Di 

3.0 A 

D 2 Di 

4000 

49.23 

6.10 

2.49 

4000 

78.9 

78.1 

78.5 

75.6 

72.1 

66.1 

4500 

29.73 

4.34 

2.51 

4500 

93.3 

94.1 

94.0 

88.6 

77.2 

57.1 

5000 

6.11 

1.08 

1.02 

5000 

72.6 

66.8 

54.3 

37.1 

25.0 

13.1 

5500 

2.86 

0.55 

0.56 

5500 

49.8 

36.2 

23.9 

13.9 

9.9 

5.0 

6000 

1.56 

0.29 

0.35 

6000 

30.4 

19.4 

11.6 

7.1 

5.3 

3.0 

6500 

0.88 

0.14 

0.23 

6500 

16.0 

10.5 

5.7 

3.9 

2.7 

1.6 

7000 

0.57 

0.08 

0.15 

7000 

9.0 

7.2 

3.3 

2.9 

1.6 

1.2 


summarized in Table Clearly, there is an effect related to the 
CLV, which is stronger for lower effective temperatures for the 
specified bands. We attribute this trend to the increasing width of 
the Fraunhofer Nat lines, when effective temperature decreases. 


ICharbonneau et aT] ( |2002| l carried out the analysis for 
HD 209458 b, which orbits a host star with an effective tem- 
perature of 6065 K (|Torres et al. 2008)1. The values given by 


Charbonneau et al. 


__ ( |2002| l read T52 x 10-^ 0.39 x 10 ^ and 

0.47 X 10“^ for the narrow, medium, and wide bands. These val¬ 
ues are compatible with the values we derived for a host star with 
an effective temperature of 6000 K. We attribute the small devi¬ 
ation to the difference in the system geometry and the different 
model spectra. 


4.2.2. Difference curves with narrow feature intervals 


From our simulations, we obtained the feature light curve, //(f,), 
by integrating the spectral signal in bands with half-widths of 
0.1 A, 0.375 A, 0.75 A, and 1.5 A, centered on either of the two 
sodium lines. The reference light curve is defined as the median 
flux obtained across a broad (5820 - 5950 A) wavelength range. 
The resulting DCs are shown in Fig. ^for stars with effective 
temperatures of 4000 K, 5000 K, and 6000 K. The associated 
DCEs are listed in Table |2l For all stars and bands considered 
here, we obtain positive DCEs. Therefore, the sodium lines show 
net brightening during transit. 

The actual choice of a reference interval in the simulated 
spectral region remains quite insignificant for the resulting DCs. 
Although there is some differential limb darkening even in the 
continuum, the level of change is much smaller than the devia¬ 
tions seen in the spectral lines. Eor instance, the DC of the two 
potential reference intervals 5860 - 5865 A and 5920 - 5940 A 
on the blue and red side of the sodium lines shows variation in 
the +2 X 10“^ range for a star with Tefr = 4000 K, which is two 
orders of magnitude smaller than that observed in the Na i lines. 

The DCs involving the sodium line show variability on a 
level of up to 10“^. The behavior is typical for transit light curves 
obtained from a limb-brightened feature (cf. Sect. |3.3| ). The de¬ 
tails obviously depend on the choice of the extraction bands and 
CLV of the stellar lines. Nonetheless, the amplitude of the curves 
shown in Eig. can clearly be sufficient to be relevant for the 
study of exoplanetary atmospheres, whose signals are expected 
to be on the order of 10“^ - 10 (e.g., Seager & Sasselov 2000| l. 


4.2.3. Treatment of the CLV-induced in-transit excess 

Through the partial occultation of the stellar disk during transit, 
the stellar CLV is imprinted on the DC. In the cases considered 
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Table 3. Properties of HD 189733. 


Property 

Value 

Source" 

Teff [K] 

The host star 

5040 + 50 

T 

log(g[cms- 2 ]) 

A cq'7+0.014 

T 

[Ee/H] 

-0.04 + 0.08 

B/T 

ysin (0 [kms^'] 

3.5 + 1 

B 

M, [Mo] 

0.806 ± 0.048 

T 

Rs [Ro] 

0.756 + 0.018 

T 

The planet and its orbit 


SMA* [Rs] 

8.81 + 0.06 

T 

Rp/Rj 

0.15463+0.00022 

T 

^orbit [deg] 

85.58 + 0.06 

T 

To [BJDurc] 

2454279.436714 

A 

period [d] 

2.21857567 

A 


Notes. ° B |Bouchy et al.|p005 i, T flbrres et al.^^2008^ , A |Agol et al.| 
1 2010 1 ; * Orbital semi-major axis 


here, it causes apparent net emission in the feature bands cover¬ 
ing the Na I Di and D 2 lines. This has to be taken into account in 
the measurement of in-transit excess absorption due to the plan¬ 


etary atmosphere as has been shown, e.g., by Charbonneau et al. 
( 2002| l and |Sing et al. (2008 1 . In particular, the corresponding 
value should be added to the measured transit depth, if we have 
not previously taken the CLV into account. 


As the DC is time dependent, the result of the averaging re¬ 
quired to obtain the DCE (see Sect. |3.2| i depends on the sam¬ 
pling of the transit. An accurate correction of the CLV-induced 
effect is only possible, if all boundary conditions like the tempo¬ 
ral sampling, instrumental resolution, and sensitivity are taken 
into account in the modeling. 


5. Observation of the center-to-limb variation in 
HD 189733 

In 2012, we obtained high-resolution transit spectroscopy of 
HD 189733 with UVES. The properties of HD 189733 are sum¬ 
marized in Table We generated synthetic spectra of the star 
based on a Kurucz model atmosphere with an effective temper¬ 
ature of 5000 K, a surface gravity of log( 0 ) = 4.5, and solar 
abundances. Again, we used our code to simulate a spectral time 
series based on the orbital configuration of HD 189733 system 
from which we derive model light curves to be compared with 
the observations. 
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Fig. 5. Simulated DCs for the Nat lines (D 2 in the left column and D| in the right column) for effective temperatures of 4000 K, 5000 K, and 
6000 K. Normalized feature light curves where extracted in four bands with half-widths 0.1 A, 0.375 A, 0.75 A, and 1.5 A centered the respective 
sodium line. The resulting DCs are plotted using solid (black), dashed (red), dotted (green), and dash-dotted (blue) lines in the same order. 


5.1. Observations and data reduction 

On July 1, 2012, we observed HD 189733 with UVES mounted 
at the VLT-UT2 (Kueyen). Between 04:06 and 08:40 UT, we 
obtained 244 spectra, whereof the first 29 were exposed for 30 s 
and the remaining for 45 s. We applied the dichroic beam splitter 
(dic2) in the 437 -H 760 nm setup. The slit width was 1" for the 
blue and 0.7" for the red arm. To optimize the temporal cadence, 
we used the “ultra fast readout” option. 

To reduce the echelle spectra, we used the UVES-pipeline in 
version 5.2.0. The pipeline performs bias subtraction, hot-pixel 
rejection, order definition, order extraction, flat holding, and the 
wavelength calibration using ThAr frames. Einally, the echelle 
orders are merged into a single Id-spectrum and Hux calibrated. 
In this work, we use the spectra from the blue and lower red chip, 
which were obtained using optimal extraction. 

We obtained reference spectra for Hux calibration during our 
observing night, but found that the spectra calibrated using the 
master response Hle^show a lower degree of large-scale struc¬ 
ture not attributable to the stellar spectrum, while the small-scale 
noise properties of the resulting spectra are identical. The Hux 
obtained with the master response files was generally lower by 
about 15 - 20% on the lower red chip, which gives an idea of 


^ see http://www.eso.org/observing/dfo/quality/UVES/qc/ 
SysEffic_qcl.html We used UV_MRSP_041130_BLUE437 for the 
blue, UV_GMRE_070511A_inaster_response_REDL76S for the lower 
red, and UV_GMRE_09O701A_master_response_REDU760 for the up¬ 
per red chip. 



-0.04 -0.02 0 0.02 0.04 

Orbital phase 

Fig. 6. Temporal evolution of seeing (top) and air mass (bottom). The 
vertical dashed lines indicate the duration of the transit. The temporal 
offset. To, refers to MIDhjd.utc = 56109.26164. 

the involved uncertainties. In our analysis, we prefer to use the 
spectra calibrated using the master response files. 

The blue and red chip are read out separately and, in princi¬ 
ple, asynchronously. Analyzing the start times as determined by 
the header keyword MJD-OBS, we find that, on average, the ob¬ 
servation obtained with the blue chip starts 0.83 + 1.1 s later than 
on the red chip. In the most extreme case, the observation ob¬ 
tained with the blue chip is delayed by 9 s compared to the red 
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chip. As the difference is, however, larger than 3 s in no more 
than six cases, the observations can be treated as simultaneous 
where not explicitly stated otherwise. 

The data cover the primary transit plus some pre- and post¬ 
transit time. In Fig.[^ we show the temporal evolution of the see¬ 
ing (header keyword HIERARCH ESO TEL lA FWHM) and the air 
mass during our observations. The duration of the transit is ex¬ 
plicitly indicated. The seeing is between 0.5” and 0.6” during the 
first half of the observations and further improves in the second 
half. 


5.2. Accuracy of the flux-calibration 


The flux calibration applied by the UVES pipeline takes the ex¬ 
posure time, gain, binning, atmospheric extinction, and air mas^ 
into account. To verify that reasonable fluxes are produced, we 
integrated the flux obtained with the individual chips and com¬ 
pare the numbers thus obtained with a model spectrum and a 
black body. The model spectrum is based on a Kurucz model 
atmosphere with an effective temperature of 5000 K, a log(g) 
value of 4.5, and solar metallicity (Castelli & Kurucz|2004[[Ku~ 


rucz|1970)l. The spectrum was synthesized using the spectrum 


program ( |Gray & Corbally||1994| l. To calculate the blackbody 
spectrum, we applied a temperature of 5040 K (see Table 
In Table |4] we list the minimum and maximum observed flux 
on each chip along with the flux predicted by the Kurucz atmo¬ 
sphere and the black body. 


Table 4. Measured and predicted fluxes at Earth. 


Chip Wvl. range Obs. flux K“ 

[A] [10“® ergcm^^ s“'] 


blue 

3732.1 

- 4999.7 

2.3-2.8 

2.9 

3.3 

redl 

5655.1 

-7595.1 

4.8-6.0 

5.6 

5.9 

redu 

7564.3 

- 9463.9 

3.5-4.2 

4.0 

4.4 


Notes. “Kurucz atmosphere (K), black body (BB). 


Although the spectral signal shows an apparent drop, which 
is remarkably well aligned with the transit, the depth of this hy¬ 
pothetical transit is about 5% and, therefore, it is too deep to 
be (solely) caused by the planet, which shows a transit depth of 
about 3% in the B filter ( [Bouchy et al. |20()5 1. While the spec¬ 
tral and acquisition signals can be combined to obtain a curve 
showing a transit of proper depth, the uncertainty in the scaling 
prevents a meaningful quantitative analysis. 

The observed flux and the prediction based on the Kurucz at¬ 
mosphere agree on the 20% level, which seems reasonable con¬ 
sidering the systematic uncertainty in the flux calibration (see 
Sect. 5.11. The residual variation is largely attributable to slit 


losses with some residual influence of air mass and atmospheric 
conditions that are likely also present. 


5.3. Wavelength stability 

In the red part of the spectrum (4 > 6000 A), telluric lines yield 
a radial velocity calibration with an accuracy of about 20 ms“' 


(e.g., Balthasar et al.|[T982 Caccin et al. |1985[ Gray & Brown 

2006| l. We used the telluric lines to study the stability of the 

wavelength calibration in our data. First, we obtained a nominal 
telluric transmission spectrum using the Line-By-Line Radiative 
Transfer Model (LBLRTM) code ( [Clough et al.|2005 1 . Second, 
we selected a total of six wavelength intervals between 6800 and 
9200 and third, we determined the (apparent) radial velocity 
shift of the telluric lines using a cross-correlation. 

Figure|^shows the resulting shift of the telluric oxygen lines 
between 6864 and 6968 A. The curve shows a sawtooth-like pat¬ 
tern with an overall amplitude of roughly 1 km s ', which we 
found for telluric lines across the entire red spectral arm. We 
attribute the variation to instrumental effects, mainly changes in 
the positioning of the stellar disk within the slit, which also man¬ 
ifest in a shift of the location of the echelle orders on the detector 
and changes in the asymmetry of the acquisition images. The slit 
width of 0.7” in the red arm was, indeed, larger than the seeing 
disk; the same holds for the blue arm. 


In Fig. [7] we show the measured spectral flux on the blue chip 
obtained by integrating the flux-calibrated spectrum and the sig¬ 
nal measured in the blue-channel acquisition images as a func¬ 
tion of time. The latter clearly represents light having missed the 
slit. Acquisition and spectral signal are strongly anticorrelated. 
The evolution of the spectral and acquisition signal after the first 
half of the observation is clearly related to the behavior of the 
seeing (Fig. with better seeing resulting in more stellar light 
entering the slit. 


0.6 



-0.04 -0.02 0 0.02 0.04 

Orbital phase 


„ 2.9 

w 

cnT 2.8 

I 2.7 
t 2.6 

-o 2.5 
2.4 
^ 2.3 

LL „ - 


Fig. 7. Temporal evolution of the measured spectral flux on the blue chip 
(squares) and (scaled) signal in the acquisition images from the blue 
camera (crosses). Dashed vertical lines indicate the transit duration. 


^ “UVES Pipeline User Manual” Issue 22.3, Sect. 11.1.21 
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Fig. 8. Radial velocity shift of telluric lines between 6864 and 6968 A. 

To coiTcct for the instrumental shifts, taking its wavelength 
dependence into account, we used the telluric corrections deter¬ 
mined in the six individual bands and interpolated to correct at 
intermediate wavelengths. 

5 . 4 . The Rossiter-McLaughlin effect 

The Rossiter-McLaughlin effect is an apparent radial veloc¬ 
ity shift caused by a change in the spectral line profile during 
the planetary transit (e.g., |Rossiter||1924| |McLaughlin||1924) l. 

^ We used the intervals 6864-6968 A, 7160-7221 A, 7221-7351 A, 
8128 - 8180 A, 8949 - 9037 A, and 9055 - 9138 A. 
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The Rossiter-McLaughlin curve provides important information 
about the planetary system, most notably, the relative orientation 
of the (sky-projected) stellar spin axis and the planetary orbit 
normal (e.g., Ohta et al.|2005[ ). 

In HD 189733, the amplitude of the Rossiter-McLaughlin 
effect is about 40 m s“', which is a factor of 25 below the instru¬ 
mental velocity shifts seen in the telluric lines. To measure the 
Rossiter-McLaughlin effect, we determined the stellar radial ve¬ 
locity by cross-correlation with a synthetic spectral model in the 
6620-6700 A, 6700-6840 A, 7780-7820 A, and 8700-8750 A 


bands. The stellar radial velocity was then corrected by subtract¬ 
ing the telluric correction determined in Sect. |5.3| and adding 
the barycentric correction. The spectral bands used in the cross¬ 
correlation were selected because they are close to the telluric 
bands used to correct the instrumental radial velocity shift, but 
are themselves only weakly affected by telluric absorption. The 
resulting curve is shown in Fig. [9] along with a nominal model 
calculated using the parameters from |Triaud et al.| ( |2009| l; [0hta| 
|et al.| ( |2005[ ). The errors were estimated from the scatter of the 
radial velocity measurements obtained in the individual spectral 
bands. Note that we did not carry out a fit. 



Orbital phase 

Fig. 9. Rossiter-McLaughlin radial velocity curve along with nominal 
model with parameters from |Triaud et al.| p009| l (top) and residuals 
(bottom). 


The Rossiter-McLaughlin effect is clearly seen. Toward the 
end of the observation run, there seems to be a systematic offset 
between model and measurement, which may due to a short¬ 
coming of the instrumental correction. Otherwise, the Rossiter- 
McLaughlin effect is compatible with the curve reported by |Tri-| 
|aud et al.1 ( |2009| l. During the transit, we obtain correlated residu¬ 
als not unlike those seen in previous works (see|Winn et al.|2006t 


Triaud et al. 


Triaud et al. 


2009|l. The residuals seen by Winn et aTj ( 2006| l and 
(|2009|l show a remarkable point symmetry with re¬ 


spect to mid-transit time, which is also indicated by our residu¬ 
als. 


5.5. Line profiles and instrumental resolution 

Changes in the line profile potentially interfere with all analy¬ 
ses based on individual line characteristics, in particular, the line 
depth or flux observed around the center of the line. To study 
the temporal stability of the line profile, we selected a total of 
19 sufficiently isolated, telluric spectral lines between 5885 A 
and 7265 A. All the lines were modeled using a Gaussian pro¬ 
file. Assuming that the lines are intrinsically narrow, the instru¬ 


mental resolution is obtained by the inverse of the FWHM times 
the wavelength. Figure show the equivalent-width weighted 
mean instrumental resolution. The temporal evolution of the res¬ 
olution is anticorrelated with the seeing (see Fig.|^. 



Fig. 10. Equivalent-width weighted instrumental resolution obtained 
from 19 telluric lines between 5885 and 7265 A. 

According to our analysis, the mean instrumental resolution 
is about 60000, but it clearly varies during the observing run. 
The behavior can likely be attributed to a change in slit illumi¬ 
nation, with the seeing disk moving in the slit. 


6. Center-to-limb variation observed in the Can H 
and K iines 

The Can H and K lines at 3933.66 A and 3968.47 A are the 
widest and strongest spectral lines in the wavelength range cov¬ 
ered by our UVES observations. These lines are expected to 
show a distinct pattern of CLV across their profiles, and their 
width makes them favorable targets for studying the CLV. Prior 
to any further analysis, we cross-correlated the observed spec¬ 
tra with a model spectrum to correct for the instrumental radial- 
velocity shift (see Sect. 5.3 1 . In the following, we adopt the 
reference continuum bands Ch (3891. 67 - 3911.67 Al a nd Ck 
(3991.067 - 4011.067 A) also used by|Melo et al.]|2006 1 . 


6.1. The line core 

The cores of the Ca ii H and K lines are highly sensitive to chro¬ 
mospheric emission (e.g., |Baliunas et al.||1995[ ). To obtain the 
DC, we integrated the spectral signal in 1 A wide bands centered 
on the Can H and K lines cores and derived two normalized 
feature light curves LCKih) and LCniti). The integration of the 
signal in the two aforementioned reference bands yields two nor¬ 
malized reference light curves, LCc^di) and LCcudi)- Based on 
these light curves, we obtained the DC as 


_ {LCKdi) + LCndi)) (LCcKdi) + LCci,(_ti)) 

^^^COre\U) — rx rx -> 


where the first term is the mean feature light curve and the sec¬ 
ond the mean reference light curve. In Fig. 11 we show the re¬ 
sulting DC, which shows a complex temporal pattern of variabil¬ 
ity. Most prominently, a steep rise is seen at about mid-transit 
time, which is followed by a decay phase, lasting for the rest 
of the observation. This behavior in chromospherically sensitive 
lines is typical of solar and stellar flares (e.g., |Fang et al.|1992[ 
IFuhrmeister et al.|2008| l. Therefore, we attribute it to a flare. Su¬ 
perimposed on the decay phase of the flare is another smaller, yet 
otherwise similar, event. We interpret this as another, potentially 
related, flare. Indeed, multiple temporally overlapping flares are 
frequently observed on active stars (e.g.,|Lalitha et al.|2013|l. In 
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Fig. 11. Top: the observed DC extracted from the Can H and K line 
cores. Bottom: close-up of the first half of the observed DC. The solid 
line shows a photospheric model DC and the dashed line a chromo¬ 
spheric model DC, assuming homogeneous distribution of the Can H 
and K core emission. 


the following analysis, we exclude the second half of the DC 
from our modeling of the CLV. 

To model the observed behavior in the first half of the ob¬ 
servation, we derived two synthetic DCs. First, we obtained a 
photospheric model, representing the DC that would be obtained 
in the absence of any chromospheric emission from our synthetic 
spectral time series. By fitting a transit model with quadratic 
limb darkening ( |Mandel & Agol|[2002[ ) to the synthetic mean 
feature and reference light curves, we obtained limb-darkening 
coefficients of a/ = 0.27 and bf = 0.33 for the feature and 
ar - 1.08 and = -0.24 for the reference. The latter are com¬ 
patible with the values of a = 1.06 and b - -0.2 derived by 
|Claret| ( |2004[ ) for stars like HD 189733 in the Sloan w'-band. Sec¬ 
ond, we constructed a DC assuming a homogeneous distribution 
of chromospheric emission across the stellar disk, i.e., zero limb 
darkening of the feature. The reference light curve was left un¬ 
changed. To model the observed DC, we additionally added a 
second second-order polynomial to the DCs to take the global 
pattern of variability into account. The actual spatial distribution 
of chromospheric emission is probably more inhomogeneous as 
observed on the Sun ( [Llama & Shkolnik||2015| ), but we refrain 
from adopting a more complex model here. 

In the bottom panel of Fig. [m we show the first half of the 
data along with the our two best-fit photospheric and chromo¬ 
spheric models. In the fitting, only the coefficients of the poly¬ 
nomial were varied. The behavior of the observed DC during 
ingress is reminiscent of the expected behavior. In particular, an 
initial decrease followed by an increase toward the mid-transit 
phase is clearly visible. The behavior is better reproduced by 
the chromospheric model, assuming homogeneously distributed 
emission in the Can H and K emission line cores. While this is 
plausible, there is also variability on a similar scale as observed 
prior to ingress, which is, therefore, most probably unrelated to 
the transit. Although the model is plausible, it remains difficult, 
if not impossible, to distinguish between intrinsic variability and 
CLV-induced changes in the cores of the Can H and K lines. 

6.2. The line wings 

In contrast to the line cores, the wings of the Can H and K 
lines are not expected to be heavily affected by stellar activity. 
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Fig. 12. Combined DC of the Can H and K lines wings along with 
our hest-fit second order polynomial model (dashed) and our best-fit 
model including the synthetic DC. The bottom panel gives the residuals 
with respect to our DC model. The solid (red) line shows the residuals 
smoothed hy a running mean with a width of eight data points. The 
vertical lines indicate the transit duration and the three top hars indicate 
the phase intervals used in the analysis of the spectra. 


In Fig. 12 we show the DC derived from the Can H and K lines 
of HD 189733. Here, our feature light curve was extracted by 
integrating the spectrum in a region between 3 A and 5 A from 
either line core, while the cores themselves were excluded. A 
separation of 3 A from the Ca ii H core also excludes the He line, 
which is not unexpectedly prominently seen during the flare (see 
Fig. [T3 ] i. We obtained four normalized feature light curves cor¬ 
responding to the 2 A wide bands on the blue and red side of 
the Ca II H and K line cores and obtained the mean feature light 
curve by averaging (cf. Eq.0. 

The DC in the line wings shows the general pattern of vari¬ 
ability expected from a limb-brightened feature light curve (see 
Sect. |3.3| l. However, there clearly is additional large-scale evolu¬ 
tion in the DC, which is reminiscent of the behavior of air mass 
(see Sect. HD- After having estimated the noise from the scatter 
in the curve, we modeled the curve, first, using a second-order 
polynomial and, second, using the polynomial plus a synthetic 
DC. The resulting models are also shown in Fig. For the 
polynomial alone, we obtained a reduced value of 1.31. Af¬ 
ter adding the model DC, we obtained a reduced value of 
1.02. Because we have not added a free parameter, the second 
model including the synthetic DC is clearly superior. The prob¬ 
ability of obtaining a reduced x^ value of 1.02 or higher from 
random noise is 40 %, which underlines the formal plausibility 
of the presented model. While there may be some evolution not 
accounted for by our model in the residuals shown in Fig. 12 in 


particular, during the first half of the transit, we consider them 
acceptable, given that no DC parameters were actually fitted 
(only the underlying polynomial). Finally, we checked that com¬ 
patible results can be produced using various feature and refer¬ 
ence bands to exclude a chance finding based on our particular 
choice of the integration bands. 


6.3. Difference spectra around Ca ii H and K 

In our next step, we aim to trace back the pronounced signal in 
the DC to changes directly observable in the spectra. Therefore, 
we selected the spectra showing the largest differences as indi- 


17 































































S. Czesla, T. Klocova, S. Khalafinejad, U. Wolter, J.H.M.M. Schmitt: The CLV across the Fraunhofer lines of HD 189733 



Wavelength [A] 

Fig. 13. Close-up of the difference between the coadded spectra in the 
periods r 2 and c. The He line is clearly seen. 


cated by the DCs. In particular, we defined the three orbital phase 
intervals ri = -0.011 + Ap, c - 0.0 + Ap, and ra = 0.011 + Ap 
with Ap = 0.0027, which results in a length corresponding to 
16% of the transit duration. The interval c covers the transit cen¬ 
ter and ri 2 cover the DC minimum before and after mid-transit 


time; they are indicated in Fig. 12 


We then coadded all spectra matching the phase require¬ 
ments, taking their relative radial-velocity shift into account, 
and then normalized the result by the average of the two mean 
flux values obtained in the Ch and Ck reference intervals. We 
then subtracted the coadded mid-transit spectrum from either of 
the reference spectra. Finally, we chose similarly phased spec¬ 
tra from our time series of synthetic spectra and carried out the 
same analysis there to obtain a model. 

Figure 13 shows a 5 A wide excerpt of the difference spec¬ 
trum obtained by subtracting the mean spectrum in interval c 
from the r 2 spectrum. As a result of the flare, a clear excess is 
seen in both the core of the Ca ii H and the He line at 3970.08 A. 


.3 -3 




Fig. 14. Top: difference of the normalized spectra obtained in time inter¬ 
vals ri and c (squares). The solid line shows the sum of our best-fit sinu¬ 
soidal plus linear model components. Vertical dashed (green) lines indi¬ 
cate the start and vertical dotted (red) lines indicate the end of an echelle 
order. Bottom: the residuals obtained after subtracting the model shown 
in the top panel from the data. The solid (red) line shows a smoothed 
synthetic difference spectrum. 

On larger scales, the difference spectrum shows structure 
that is not attributable to the star. In the top panel of Fig. 
we show the difference spectrum resulting from the comparison 
of the coadded precenter reference spectrum and central spec¬ 
trum. The difference spectra were rebinned on a 0.5 A grid to 
improve visibility. Clearly, the difference spectrum shows a pe¬ 


Table 5. Best-fit parameters for fit of difference spectrum 


Parameter 

Best-fit values 


ri - c 

r 2 - c 

A 

1.6 X 10-^ 

2.5 X 10-4 

</> 

1.82 

1.82 (fixed) 

xo 

1.5 X 10-4 

-2.3 X 10-4 

0 [A-'] 

3.3 X 10-® 

-9.3 X lO-'^ 


riodic, sine-like signal, whose period of 34 A is compatible 
with the spacing of the UVES echelle orders in this spectral re¬ 
gion for our instrument setup. Therefore, we argue that this sig¬ 
nal is probably related to the blaze correction and order merging 
of the spectra. 

Close inspection of the synthetic and observed difference 
spectrum shows that many details are not well reproduced. It 
must be noted that UVES resolved hundreds of spectral lines in 
the region, which all show a distinct CLV. To test whether the 
overall behavior of the observed difference spectrum matches 
the simulation nonetheless, we smoothed the synthetic difference 
spectrum using a running mean with a window size of 2.5 A and 
compared the smoothed spectrum to the observation. In partic¬ 
ular, we fitted the observed difference spectrum using a model, 
which is the sum of a sinusoidal component, a linear component, 
and the synthetic difference spectrum, i.e., 

r{A) = A sin(27rv + ^) + xo + gx(A- 4000 A) H- SYD. (10) 

Here, SYD represents the smoothed synthetic difference spec¬ 
trum. The linear and sinusoidal model components represent 
changes unrelated to the occultation of the stellar disk. After fix¬ 
ing the period, v“', of the sine to 34 A (the order spacing), we 
were left with four fit parameters: the constant, xq, the gradient, 
g, the amplitude. A, and the phase, (p, of the sine. Our best-fit 
parameters are given in Table During the fit, we excluded the 
+2 A wide region around the Ca ii H and K line cores and the 
same region centered on the He line, which is clearly seen dur¬ 
ing the flare (see Eig. [T3J. 

The solid curve in the top panel of Eig. shows the sum 
of the linear and sinusoidal components from our best-fit model 
along with the observed difference spectrum. In the bottom 
panel, we plot the residuals with respect to the curve shown in 
the top panel along with the smoothed synthetic spectrum, which 
is the last model component. While the scatter remains large, the 
overall behavior follows the expectation. 

In Eig.[^ we show the result of the same experiment using 
the post-center reference spectrum. Again, we fitted the same 
model to the residuals. As the sinusoidal component is less pro¬ 
nounced here and we argue that it is related to the order merg¬ 
ing, we fixed the phase. In this case, we obtain a stronger linear 
component (see Table |^, which we attribute to a change in the 
atmospheric transparency not completely removed by the flux 
calibration. 


6.4. Influence of the systematic effects in the difference curve 

The analysis of the difference spectra in Sect. |6.3| clearly demon¬ 
strated the presence of systematic effects in the spectra. To quan¬ 
tify the impact of these effects on the DC derived in Sect. 6.2[ we 
simulated spectra, including both a sine wave with variable am¬ 
plitude and a linear term with variable gradient. The phase of the 
sine was fixed to the value obtained from the fit to the residu¬ 
als because this is not expected to change for an effect related to 
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Wavelength [A] 


Fig. 15. As in Fig.y^for the time intervals r 2 and c. 


the echelle orders. Based on these simulations, we quantify the 
impact on the DC 


ADC 

AA ■ 


-0.04 and 


ADC 

Ag 


= -54 A . 


( 11 ) 


Given the measured amplitude of 1.5 x 10“^ of the sine wave, 
we estimate a resulting change on the order of 10 in the DC. 
The expected change in the DC resulting from the gradient is 
54 Ax 10“^ A“' = 0.54X10“^, where 10“^ A“^ is on the order of 
the observed value. The change in the gradient should, however, 
be slow because it is probably related to a change in air mass. We 
speculate that the global trend seen in the DC, shown in Fig. [T^ 
is related to this effect. 

Although there are instrumental effects, we argue that their 
magnitude remains insufficient to explain the observed DC (see 
Fig.[g. Furthermore, the temporal association with the transit 
seen in the DC makes an instrumental origin unlikely. In prin¬ 
ciple, nonlinearity of the CCD could cause an effect temporally 
associated with the transit because the star appears necessarily 
darker when the planet traverses the stellar disk. However, the 
signal seen in the DC is not proportional to the overall light 
variation expected during the transit and the signal observed on 
the CCD is on the order of 1000 ADUs in the region of the 
Can H and K lines, which is far below the saturation limit of 
65 000 ADUs. Therefore, we argue that nonlinearity of the CCD 
is unlikely to be the origin of the observed change in the DC. 


7. Center-to-limb variation observed in the Nai Di 
and D 2 lines 

Our analysis of the Can H and K lines shows flaring activity 
affecting the line cores from about mid-transit time. To study 
the effect on the Nai Dj and D 2 lines, we created mean spec¬ 
tra in three time intervals covering orbital phases from -0.02 to 
-0.01, -0.01 to 0.0, and 0.01 to 0.02. While the first two inter¬ 
vals do not cover the flare, the third interval covers its peak phase 
as observed in the Can H and K lines (cf. Fig. [TT|l. In Fig.fThl 
we plot the ratios of the averaged spectra in the three intervals. 
While the Nai Di and D 2 line cores show a clear fill-in during 
the flare, the ratio of the out-of-flare spectra shows no strong 
change in the core region. The core is probably also affected by 
the planetary atmosphere (e.g., Wyttenbach et al. |2015| l, which 
we do not consider here. To minimize the impact of the flare on 
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Wavelength [A] 

Fig. 16. The solid (red) line shows the ratio of the average peak-flare 
spectrum (orbital phases 0.01 to 0.02) and the first out-of-flare reference 
spectrum (phase -0.02 to -0.01). For comparison, the dashed (black) 
lines shows the ratio of the averaged spectra obtained during the second 
and first out-of-flare reference phases, which shows no strong change. 
The dotted (blue) line shows the scaled stellar spectrum for reference. 


our CLV analysis, we excluded the line core. In particular, we 
disregarded the central ±0.3 A region of the lines, where the line 
profile is dominated by the central Gaussian component. The re¬ 
gion is also indicated in Fig.[T6| 
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Fig. 17. Region covering the Nat Di and Do lines in our first spectrum. 
Top: shaded (red) regions show the wide reference bands and squared 
regions the narrow spectral bands. Bottom; close-up of the Nai Di and 
Do lines. The shaded (red) regions have half-width of 0.3, 0.6, 0.9, and 
1.5 A. 


In Fig. we show the observed, flux-calibrated spectrum 
around the Nai Dj and D 2 lines. Note that neither a barycen- 
tric correction nor any further normalization was applied to the 
spectra. From the observed series of spectra, we extract feature 
light curves using three different feature bands with half-widths 
of 0.6 A, 0.9 A, and 1.5 A centered on the cores of the sodium 
lines. These core regions, however, were disregarded by exclud¬ 
ing the inner ±0.3 A region. We then obtained two feature light 
curves for both line cores by individually integrating the spectral 
signal in the blue and red wing. 


In our calculations, we used both wide (20 A) and narrow 
reference bands (see Fig.ITTll. The wide reference bands are de¬ 
fined between 5840 - 5860 A and 5925.875 - 5945.875 A. To 


test the influence of the reference bands on the outcome, we de¬ 
fined two further pairs of wide reference bands by shifting both 
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bands toward the Nat Di and D 2 lines in two steps of 10 A. Ad¬ 
ditionally, we defined narrow reference bands (5868 - 5872 A 
and 5910.5 - 5912.5 A), which contain no strong spectral lines 
(Fig. 17 1 . To correct for the radial-velocity shift of the individual 
spectra, we carried out a cross-correlation with a model spectrum 
in the 5850 - 5870 A band and adjusted the start- and end wave¬ 
lengths of the integration bands according to the radial-velocity 
shifts thus derived. With this approach, we account for all present 
shifts, whether they are instrumental or physical in origin. 


The resulting DCs are shown in Fig. 18 for both line cores 


along with model curves obtained from our simulations of spec¬ 
tral time series. Around the center of the transit, we find the be¬ 


havior of the curves shown in Fig. 18 to be consistent with our 
model DCs. In particular, the DCs of both sodium lines and all 
chosen feature bandwidths show a drop during in- and egress 
phases and a peak at transit center. The synthetic curves show 
that the CLV effect is not identical for both sodium lines, but a 
little stronger for the D 2 line for the chosen integration bands. 
This reflects a difference in width and strength of the Nat Di 
and D 2 lines, with the D 2 being a little broader and deeper in 
stars like HD 189733. 


At orbital phase 0.03, the observed DCs show a drop for 
which we currently have no definite explanation. We speculate 
that this drop might be related to the shift in the curve of the 
Rossiter-McLaughlin effect seen toward the end of the observa¬ 
tion. At this point, the instrumental resolution and radial-velocity 
shift also show a general change in their behavior (see Figs.|^. 
Whether the evolution in the DC is ultimately related to an in¬ 
strumental effect or a true, physical change in the stellar spec¬ 
trum, e.g., due to the flare or a change on the surface caused by 
a spot remains unknown. 


In Fig. 19 we show the combined DC from both Nat D] 
and D 2 line cores (0.3 - 0.6 A band) up to an orbital phase of 
0.03 along with a linear model and a linear model plus our syn¬ 
thetic DC. The curves were produced using the narrow reference 
band. However, we find little change if other reference bands are 
used. The error was estimated as the scatter in the observed DC. 
Clearly, the synthetic model DC explains the general behavior of 
the observation. In particular, we obtain a value of 369 using 
only the linear model and a value of 279 if the full model is used. 
In both cases, only the two parameters of the linear model were 
fitted. The profile of the synthetic DC is completely determined 
by the transit geometry and the stellar model atmosphere. There¬ 
fore, we are left with 195 degrees of freedom. It seems, however, 
that the observed DC shows further structure not accounted for 
by the model. Before and after mid-transit time, there are two 
drops in the residuals. The structure is, however, not entirely 
symmetric with the deviation stronger during the first half of the 
transit. 


To study the effect of variable telluric absorption on our mea¬ 
surements, we used the molecfit code to correct our spectra for 
telluric contamination ( |Smette et al.|20l5||Kausch et al.|2015| l. 
The DCs produced using the spectra corrected for telluric ab¬ 
sorption do not significantly differ from those that are not cor¬ 
rected for telluric absorption. The main difference is a system¬ 
atic pattern changing in phase with air mass. The amplitude of 
this pattern is, however, only 0.5 x 10“^. Formally, we obtained a 
reduced value of 1.48 for the CLV model with telluric correc¬ 
tion, which is even a little higher than the value of 1.43 obtained 
without telluric correction. As the CLV-related wave-like pattern 
remains essentially unaffected, we conclude that variable telluric 
absorption is unlikely to be the origin of the observed behavior. 



Orbital phase 


Fig. 19. Top: DC observed in the 0.3 - 0.6 A band around both Nai Di 
and D 2 cores with respect to narrow reference intervals. The dashed 
(red) curve shows a linear fit and the solid (black) line a linear fit plus 
our synthetic DC model. Bottom: residuals with respect to the synthetic 
model. The solid (red) curve shows the residuals smoothed using a run¬ 
ning mean with a width of eight data points. 


Finally, we verified that the change in instrumental resolution 
detected in the UVES spectra (see Sect. 5.5 1 is not responsible 
for the observed behavior. In particular, we simulated DCs ap¬ 
plying instrumental resolutions of 55 000 and 65 000 represent¬ 
ing the observed change in resolution (see Sect. 5.51. The result¬ 
ing change in the DCs remains too small to be of any relevance 
in the current study. 


8. Summary and discussion 

We present DCs of the Can H and K and Nai Di and D 2 lines 
derived from transit spectroscopy of HD 189733 obtained with 
UVES. Our DCs clearly show the effect of differential CLV in 
the wings of the Can H and K and Nat Di and D 2 lines of 
HD 189733. 


8.7. The center-to-limb variation 

Our analysis of the Ca ii H and K line cores revealed strong vari¬ 
ability. Most notably, an increase in the strength of the cores, 
which we interpret as a stellar flare, was observed at about mid¬ 
transit time. We found that the core DC obtained before the flare 
is better reproduced by a model assuming no limb darkening in 
the Ca II H and K line cores than by a model with photospheric 
limb darkening. This finding is compatible with a homogeneous 
brightness of the chromospheric core emission across the stellar 
disk. Yet, the actual spatial distribution is probably more com¬ 
plex ( [Llama & Shkolnik|2015| l. The Ca ii H and K emission line 
cores also show variability before the ingress, which is proba¬ 
bly unrelated to the transit. While an earlier ingress either due to 
planetary material or an extended stellar chromosphere is con¬ 
ceivable (e.g., [Czesla et al.|2012 Haswell et al.|201^ , intrinsic 
variability, independent of the transit, appears to be a more plau¬ 
sible interpretation in view of the flare observed later. In any 
case, it remains difficult to distinguish between effects caused by 
the transiting planetary disk and intrinsic stellar variability in our 
data set. 

To limit the impact of activity, we focused our analysis on the 
wings of the Eraunhofer Can H and K and Nai Di and D 2 lines. 
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Fig. 18. Difference curves for Nai Do (left column) and Di (right column) for three feature bands with half-widths of 0.6 A, 0.9 A, and 1.5 A (top 
to bottom). Top panels show the data and the model (solid, red line), and the bottom show the residuals. The line core (±0.3 A around the line 
center) has been excluded. The reference bands were 5840 - 5860 A and 5925.875 - 5945.875 A. The effect of the CLV is different in the two 
lines and it is weaker for broader spectral bands. 


These are expected to be less strongly influenced by activity, but 
are still sensitive to the stellar CLV. The shapes of both the DCs 
based on the Ca ii H and K and Na i Dj and D 2 line wings clearly 
show the behavior expected from the stellar CLV. In particular, 
we find the line wings to be limb brightened compared to the 
reference continuum. 


The observed DCs are well reproduced by a model composed 
of a synthetic DC and a low-order polynomial. Including the syn¬ 
thetic DC model clearly improves the fit compared to a low-order 
polynomial alone. In particular, the reduced values changed 
from 1.31 to 1.02 in for the Can H and K lines and from 1.89 
to 1.43 for the Nai Di and D 2 lines. Nonetheless, we still obtain 
some correlated residuals, especially in the modeling of the DC 
derived from the Na i Di and D 2 lines. 


There may be several reasons for these residuals: first, the 
modeling may be incomplete. |Hayek et al.| ( |2012| l find better 
fits to transit light curves with CLV models based on three- 
dimensional LTE model atmospheres, but our modeling is based 
on one-dimensional LTE models. Also, non-LTE effects are ex¬ 
pected to influence the formation of the strong Eraunhofer lines 
(e.g., Lind et al.|201 l| l. Because the orbit of HD 189733 b shows 
only mino r misalignment wit h the stellar spin axis (J3 — -A — 
0 . 85 !°; 2 «°, 
feet should 


Triaud et al. 


2009 1 , we expect that this kind of an ef- 
be symmetnc with respect to mid-transit time. There¬ 
fore, atmospheric model deficiencies can hardly account for all 
the residuals. Second, there may be instrumental effects. Our 
data show a complex pattern of instrumental radial velocity shift, 
and we identified residuals, which we attributed to instrumen- 
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tal effects, in the difference spectra of the Can H and K lines. 
Clearly, instrumental effects could give rise to residuals on the 
observed scale. Between orbital phase -0.01 and 0, there seems 
to be a structure in the residuals of both the Ca ii H and K and 
Nai Di and D 2 lines, which may have a physical origin. Stellar 
activity provides viable candidates for this kind of structure. 


8.2. The influence of stellar activity on the difference curve 


Starspots and flares are prominent manifestations of stellar ac¬ 
tivity. In the following, we briefly discuss their potential impact 
on the observed DCs. Indeed, starspot occultations have been 
observed in transit light curves of HD 189733 b earlier (e.g., 
Pont et al.|2008t Sing et al.|201 11 , which makes a spot crossing 
event a plausible scenario. In particular, as we certainly observed 
HD 189733 in an active phase as evidenced by the tentative flare 
starting near mid-transit time. 

Starspots crossed by the transiting planet cause “bumps” in 
the broadband transit light curve because comparably darker re¬ 
gions of the stellar disk are occulted (e.g., |Pont et ar||2008[ 
Czesla et al. 2009 | l. Eor the DC, the height difference of fhe 
bump in fhe feature and reference light curve is relevant. As the 
wings of the Nai lines become deeper for lower effective tem¬ 
peratures (see Pig.[^, the spot appears darker in the wings of the 
Na I lines than in the surrounding continuum, given that its spec¬ 
trum can still be described by an otherwise normal photospheric 
model. This gives rise to a more pronounced bump in feature 
light curves covering the Nai line wings and, thus, a net increase 
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in the DC. The same applies to the Ca ii H and K lines. The am¬ 
plitude of the effect depends on the spectral characteristics and 
shape of the spot. For a large spot with effective temperature 
of 4000 K ( |Pont et ar]|2008| ), we expect that the spot-induced 
bump can be of the same scale as the amplitude of the actual 
DC because the change in the profile of the Nai lines is quite 
pronounced between 4000 K and 5000 K (see Fig. [^. 


Also starspots not occulted by the planet affect the DC be¬ 
cause they cause an effective wavelength-dependent change in 
the planet-to-star radius ratio and, thus, alter the normaliza¬ 


tion of the light curves, which changes the transit depth (Sing 


|et al.|2008j [Czesla et al.|2009| ). While the effect of the starspot 
crossing is local, unocculted spots modify the entire transit light 
curve, given that they do not rotate on or off the disk. At least in 
the case of HD 189733, this seems improbable, because we only 
observed a little more than one percent of the stellar rotation pe¬ 
riod. 


The flare observed as a rise in the Can H and K emission 
line cores at about mid-transit time could influence the DC both 
via chromospheric and photospheric emission. During flares, the 
strength of chromospheric emission lines including Ca ii H and K 
and Na I Di and D 2 are expected to increase (e.g., 

|Schmitt|2004| l. While we exclude the cores in the 
line wings, some contribution could still be present. A fraction of 
the chromosphere and upper photosphere is thought to be heated 
to temperatures of 10000 K ( |Haisch et al.|199i] l. Effectively, 
a part of the quiescent photospheric emission would be replaced 
by emission from a 10000 K hot spot. The spectra of hot pho¬ 
tospheres are dominated by broad Balmer lines, while lines of 
other ions (including Nai) tend to be narrow. While the exact 
impact of the flare depends on its location, size, and spectrum, 
its effect should be temporally limited, and we expect the main 
contribution to take place during the leading impulsive phase. 


Fuhrmeister & 
analysis of the 


8.3. The planetary atmosphere of HD 189733 b seen in 
Na I D] and D 2 


Difference curves have been used to detect planetary excess ab¬ 
sorption in the Nai Di and D 2 lines during transit. While there 
are residuals in our DC not accounted for by our modeling, their 
temporal behavior is rather untypical for a signal related to the 
planetary atmosphere (see Fig. [T^ . In particular, only limited 
sections before and after mid-transit time show deviations, while 
the planetary atmosphere should be visible during the entire 
transit. The presence of residuals of ultimately unknown origin 
makes it difficult to derive a reliable estimate of the atmospheric 
excess absorption in the wings of the Nai Di and D 2 lines in the 
frame of our current analysis. 


Detections of in-transit excess absorption in the Nai Dj and 
D 2 lines attributable to the planetary atmosphere have been pre- 
sented by|Redfield et al.| ( 2()08|l,|Huitson et al. ( |2012|l, and Wyf 


tenbach et al. (2015 1 . Both |Huitson et al.' ( 2012| l and Wyttenbach 

et al. ( 2015| l find the bulk of planetary atmospheric absorption 
conce ntrated in the cores of the Na i Di and D 2 lines. |Wyttenbach] 
et al. I derived a FWHM of 0.52 ± 0.08 A for their signal. As we 
exclude the inner +0.3 A region around the Nai Di and D 2 line 
cores in our analysis to minimize contamination by stellar activ¬ 
ity, we also exclude about 80 % of their signal. In the Na i Di and 
D 2 lines, they report a relative transit depth of (3.2 + 0.31) x 10“^ 
so that the maximum signal we expect is about 0.6 x 10“^. There¬ 
fore, we do not regard the lack of an obvious signal related to the 
planetary atmosphere in our current analysis to be a contradic¬ 


tion of the findings of [Wyttenbach et ar] ( |2015[ ). We also do not 
exclude a contribution by the planetary atmosphere in the DC. 

An inspection of Fig. 3 in |Wyttenbach et aT] ( 2015| l suggests 
that their curves may also show the effect of the stellar CLV. 
The curves seem to show a bump in the middle of the transit, 
which resembles the DCs we obtain for limb-brightened fea¬ 
tures. In fact, our simulation for their 12 A wide band indicate 
a peak-to-peak amplitude of about 0.06% for the DC, which 
seems to match the distribution of their residuals. In the smaller 
0.75 A band, the expected CLV-induced variation is on the or¬ 
der of 0.4 % (see Fig. El, which is again compatible with their 
findings. We emphasize, however, that their detection of excess 
sodium absorption would not be challenged by this. From our 
models, we obtain a net CLV-induced transit depth of 0.006 % 
in the 12 A wide band, which is on the order of their error. As 
this shift would be systematic, the level of sodium absorption 
derived by [Wyttenbach et al.| ( [2015| l may be underestimated by 
this rather small amount. 


8.4. The impact of the temporal sampling of the transit 

The spectrum of the transited stellar disk evolves continuously as 
the planet proceeds along its path and therefore, the DC is time 
dependent. This implies that the temporal sampling of the transit 
is highly relevant in the analysis of transit spectroscopy. If, for 
instance, only the central part of the transit of HD 189733 b was 
covered, the wings of the Na i Dj and D 2 and Ca ii H and K lines 
would actually show net emission. The reverse is true, if only the 
in- and egress phases were covered. This effect occurs without 
requiring any change in the physical radius of the planet. 


8.5. Relevance of the CLV in large planets 

If the planetary atmosphere can be treated as a thin shell en¬ 
shrouding the planet, it extends the planetary disk seen during 
transit by a small amount, ARp. The area, Aa, of the resulting 
atmospheric annulus reads 



7:(Rp + ARp)^ - nR^ 

ttR] 



- 2pAp . 


( 12 ) 


When the strength of the atmospheric absorption is proportional 
to the atmospheric area, the effect observable in the spectrum 
is proportional to the planet-to-star radius ratio, p, and the ex¬ 
tent of its atmosphere. Therefore, large planet-to-star radius ra¬ 
tios and large atmospheric scale heights are favorable to detect 
atmospheric absorption. 

A simplified model of the transit light curve with quadratic 
limb darkening can be obtained for small planets, for which the 
normalized transit light curve as a function of the limb angle, p, 
reads 


n,p(a,b,p,p) 


l-p 


2 


1 -ajl -p)-b(l -pf 
^36 / 


(13) 


Here, the limb darkening is parameterized according to Eq. 
(cf. [Muller et al.|[2013| Eq. B.l). Accordingly, the DC can be 
approximated by 


DCsp(p,p) = n,p{af,bf,p,p) - n,p(ar,br,p,p) ~ p^ , (14) 


which depends on the properties of the normalized feature (f) 
and reference (r) light curves. While this simplified model does 
not reproduce the in- and egress phases appropriately, Eq.[l4]can 
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be used to estimate the val ue of the DC in the transit center. For 
the cases discussed in Sect. 33 we obtain values of -8.9 x 10“"^ 
for the limb-darkened feature and -1-10.4 x 10“^ for the limb- 
brightened feature, which both agree well with the results ob¬ 
tained from the more elaborate models (cf. Fig.[^. 

Equation 14 demonstrates that the amplitude of the CLV- 
induced effect on the DC is proportional to the squared planet-to- 
star radius ratio, while the strength of the signal induced by the 
planetary atmosphere grows linearly with the radius ratio. Addi¬ 
tionally, smaller stars, which potentially provide more favorable 
(i.e., larger) radius ratios also tend to be cooler. In cool stars, 
however, the CLV-induced variation is also more pronounced. 
This underlines the importance of an accurate knowledge and 
modeling of the stellar CLV in the study of planetary atmo¬ 
spheres. 


9. Conclusion 


We clearly detect the CLV-induced effect in a spectral transit 
time series of HD 189733, and our findings are generally com¬ 
patible with predictions from ID-LTE models. Remaining devi¬ 
ations may be explained by more elaborate modeling of the stel¬ 
lar CLV, stellar activity, instrumental effects, or a combination 
thereof. The CLV-induced effect is stronger for cooler stars and 
its relative importance increases for larger planet-to-star radius 
ratios. Our results demonstrate that the variation in the stellar 
CLV across individual spectral lines can be detected with transit 
spectroscopy and that the effect should be taken into account 
in the study of planetary atmospheres via transmission spec¬ 
troscopy. 

Transiting planets provide a valuable opportunity to resolve 
the stellar disk, which has been used to study the surface dis¬ 
tribution of starspots (e.g., Pont et aL]|2008[ Huber et al.||2010 


|Sanchis-Ojeda & Winn|[2()l 1 1 . Here we show that the mov¬ 

ing planetary disk can also be used to study the spatially re¬ 
solved stellar spectrum, which is required to constrain three- 
dimensional hydrodynamic models of the stellar photosphere 
( |Dravins et al.|20T4l l. 

Given appropriate modeling, the prominent CLV in the stel¬ 
lar Lraunhofer lines could also be used to reconstruct the plane¬ 
tary path across the stellar disk in terms of the temporal change 
in the occulted limb angle. In its use of the spectrum, this ap¬ 
proach resembles the Rossiter-McLaughlin effect, which is fre¬ 
quently used to study the orbit geometry of exoplanets (e.g., 
|Winn et al.|2006[|Triaud et al.|2009| l. While the strength of this 
effect depends on the properties of the stellar atmosphere, stellar 
rotation is not required to see it. Clearly, the CLV in the specific 
ensemble of stellar spectral lines used to derive the radial veloc¬ 
ity shifts should also be taken into account in the modeling of 
the Rossiter-McLaughlin effect. We speculate that the correlated 
residuals in the Rossiter-McLaughlin curve obtained by |Winn| 
|et al.| ( [2006| ), [Triaud et al.| ( |2009| l, and in this work could be at¬ 
tributable to the stellar CLV. 

As noted by |Yan et ah] ( |2015| ), the CLV-induced effect on 
the light curves derived from transit spectroscopy becomes more 
pronounced when narrow spectral bands are used for integration 
in strong spectral lines. Although the use of broader bands re¬ 
duces the influence of the CLV, a number of studies have found 
the planetary atmospheric signal to be concentrated in the cores 
of the line. Therefore, the use of broader bands may not always 
be desirable. At any rate, an accurate treatment of the stellar CLV 
is indispensable in the analysis of transit spectroscopy of exo¬ 
planets. 
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